Testing the Dark Matter Interpretation of the PAMELA Excess through 
Measurements of the Galactic Diffuse Emission 
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We propose to test the dark matter (DM) interpretation of the positron excess observed by the 
PAMELA cosmic-ray (CR) detector through the identification of a Galactic diffuse gamma-ray 
component associated to DM-induced prompt and radiative emission. The goal is to present an 
analysis based on minimal sets of assumptions and extrapolations with respect to locally testable or 
measurable quantities. We discuss the differences between the spatial and spectral features for the 
DM-induced components (with an extended, possibly spherical, source function) and those for the 
standard CR contribution (with sources confined within the stellar disc), and propose to focus on 
intermediate and large latitudes. We address the dependence of the signal to background ratio on the 
model adopted to describe the propagation of charged CRs in the Galaxy, and find that, in general, 
the DM-induced signal can be detected by the Fermi Gamma-ray Space Telescope at energies above 
100 GeV. An observational result in agreement with the prediction from standard CR components 
only, would imply very strong constraints on the DM interpretation of the PAMELA excess. On the 
other hand, if an excess in the diffuse emission above 100 GeV is identified, the angular profile for such 
emission would allow for a clean disentanglement between the DM interpretation and astrophysical 
explanations proposed for the PAMELA excess. We also compare to the radiative diffuse emission 
at lower frequencies, sketching in particular the detection prospects at infrared frequencies with the 
Planck satellite. 
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I. INTRODUCTION 

> n 

! Recent measurements of the positron fraction in cosmic rays (CRs) up to 100 GeV by the PAMELA experiment [Ij 

. (and of the sum of the electron plus positron fluxes by ATIC and PPB-BETS Q) have triggered a lot of interest 
on the possibility that there may be a dominant contribution to these terms from the annihilation or decay of dark 
, matter (DM) particles, assumed to compose the dark halo of the Milky Way. Indeed, the sharp raise in the positron 
' fraction detected by PAMELA above 10 GeV, confirming and super-exceeding previous hints of an anomaly 0, |^, @| 
is a feature which cannot be accommodated within the " standard picture" , with primary electrons accelerated in 
supernova remnants (SNRs) and secondary positrons produced mainly from the interaction of primary cosmic-rays 
with the interstellar medium (ISM) during propagation. It is instead suggestive of an extra primary source of 
positrons. Pulsars are well-motivated candidates for this role, see, e.g., 0,11) 3- Another interesting possibility 
is the recent proposal that positrons (and electrons) are secondary products of hadronic interaction within the cosmic 
ray sources [lOj. 

'c^ A contribution to the positron flux from dark matter was predicted long ago, see, e.g. [ill, El, [3; depending 

on the dark matter particle mass and on the shape of the electron/positron yield from annihilation or decay, the 
excess can be accommodated, see, e.g., p^ . A few puzzling issues arise however in connection to the dark matter 
interpretations: the required positron injection rate is fairly large, and, supposing the signal is due to annihilating 
particles, this points to pair annihilation rates much larger than the typical rates for thermal relic weakly interacting 
massive particles (WIMPs) and/or to substructure enhancements which should be unrealistically large. Also, it is to 
some extent surprising that a dark matter signal emerges clearly in the positron flux and it is still hiding in other 
indirect detection channels: for typical particle physics dark matter candidates (where the "typical" here stands for the 
absence of a mechanism forcing the coupling of the dark matter particle to light leptons only) signals to background 
ratios tend to be larger for the antiproton flux (and even more, in perspective, for the antideuteron flux), and for the 
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gamma-ray component produced via neutral pion decays. These features have guided recent work on both proposing 
new dark matter candidates, as well as setting constraints on the range of those that had been previously proposed. 

The issue of whether the standard cosmic ray picture could give a fair description of the electron/positron popu- 
lations in the Galaxy and of whether the DM component could give a sizable contribution had actually been risen 
a few times even before the recent refined measurements of their local fluxes. Indeed, although there is a wide 
consensus in associating the observed multi-wavelength diffuse emission in the Galaxy to CR interactions with the 
ISM and magnetic field, CR spectra and propagation parameters tuned to describe the local picture can lead to an 
underestimation of the flux in some regions of the sky. High energy electrons emit radio to microwave photons via 
synchrotron radiation on the Galactic magnetic flelds and an excess in the innermost region of the Galaxy has been 
claimed based on WMAP data (th e so-called " WMAP haze" ) when scaling the measured synchrotron radio intensities 
to the microwave band [H, [l^, Ell • Equally important is the Inverse Compton (IC) radiation component of high 
energy electrons on the starlight and microwave photon backgrounds. One of the models put forward to address the 
excess above 1 GeV reported by EGRET in the observations of the diffuse gamma-ray emission of the Galaxy [l^ 
considered the possibility of a dominant contribution from IC on starlight of an electron population in the Galaxy, on 
average, much larger than what is locally measured 20] ("optimized" Galprop model). Recent preliminary data [2lj | 
from the FERMI gamma-ray telescope (formerly GLAST), however, do not seem to confirm the EGRET excess, with 
the level of the diffuse flux at high energy {E < 10 GeV) possibly consistent with a standard cosmic-ray setup and 
pointing to an instrumental effect in EGRET. Other excesses, which are not commonly connected to a DM interpre- 
tation, has been reported in the X-ray emission at the Galactic Ridge 22] and in the TeV 7-ray emission measured 
by Milagro [23j . The first relies on a quite complicate region, where simplified assumptions on the primary injection 
spectra and on the ISM properties can break down. The contribution of unresolved source in the Milagro data can be 
large [13] . Moreover, we are mostly interested in a fair description of the local region, which involve only the diffuse 
emission at intermediate and high latitudes and longitudes. 

Most attempts to insert the dark matter interpretation of the PAMELA data into a more global picture have 
involved either other species such as gamma-ray and antiproton yields, or other dark matter environments (e.g., the 
central region of the Galaxy, or Galactic satellites, instead of the local dark matter population which, in case the 
dark matter interpretation holds, would be most likely responsible for the measured positron flux). These kind of 
comparisons are inevitably model dependent. Comparing yields is certainly very powerful when having a reference 
particle physics model in mind; on the other hand, as already mentioned, one may just tune the dark matter candidate 
in such way that only light leptons are produced as detectable yields. Radiative components are instead unavoidably 
associated to electron/positron yields: IC emission of a 100 GeV to 1 TeV electron on 1 /xm starlight photons gives 
gamma-rays with energies peaked in about the range 50 GeV to 5 TeV; the associated synchrotron emission on a 1 ^G 
magnetic field is peaked between 50 to 5000 GHz (scaling linearly with the magnetic field). Having normalized the 
electron/positron yield to the locally observed flux, the extrapolation for the radiative emission in the local portion of 
the Galaxy and its neighborhood is fairly solid, since although it is true that there are uncertainties in the astrophysical 
model involved, such as the parameters in the cosmic-ray propagation model and the level of the stellar radiation 
and magnetic fields, these need to be in turn readjusted to the local measurements (regarding the local electron flux 
as well as from ratios of secondaries to primaries in cosmic rays). Extrapolations to, e.g., the central region of the 
Galaxy to compare, e.g., with the WMAP haze, are much more uncertain since they rely on the extrapolations for the 
astrophysical parameters (with various kinds of degeneracies which can be involved) as well as for the level of dark 
matter yield (in connection, e.g., with the dark matter density profile). 

We will study the dark matter contribution to Galactic radiation seed at intermediate to large latitudes. The 
emission in such portion of the sky can actually probe the local environment and this test of the DM interpretation 
of the PAMELA excess can be thus considered self-consistent. 

We will analyze the dependence of the predicted signal on the uncertainties in the cosmic-ray propagation model. 
We will also sketch the differences with respect to the picture in which the source (s) accounting for the PAMELA 
positron fraction is confined to the stellar disc rather than being spread to on a much larger vertical extension such as 
in case of the dark matter source. The reference to test experimentally our proposal will be the FERMI gamma-ray 
telescope (and, to a minor extent, the PLANCK satellite [25]). 

Recent analyses of the diffuse Galactic emission induced by DM models in light of the PAMELA and ATIC data 
includes Refs. [H, [13, [H, [H, [s^, [IH . 

The paper is organized as follows. In Section [TTl we present the description of the particle propagation in the 
Galaxy. In Section Hill we discuss the spectral and spatial distribution of the CR and DM sources. In Section |TVl we 
present the results for the final /e~ distribution and for the multi-wavelength diffuse emission induced by CRs and 
DM. Section IVl concludes. 
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TABLE I: Benchmark models of propagation. The spectral index break for protons and electrons is at 9 GeV in the cases with 
Kolmogorov diffusion, and at 40 and 9 GeV, respectively, in the Kraichnan case. The scale of diffusion in the model B5 is taken 
to be 2s = 4 kpc. 



II. COSMIC-RAY PROPAGATION IN THE GALAXY 



We adopt the description of cosmic-rays as particles propagating in a determinate environment (i.e., disregarding 
the effects induced on the ISM by the interaction with CRs). The CR propagation equation for a particle species i 
can be written in the form [3211: 



= ^ ■ [Dxx^ni ~ VcUi) + —p Dpp- 2"'»~fl~ t>ni - -{V ■ Vc)ni +q{f,p,t)-\ -H (1) 



dt ^ ' ^ ' ' " " dp"^ dp p2 dp 



P 



3 



Tr 



where rii is the number density per particle momentum [ni{p)dp — Ni{E)dE, with Ni{E)dE being the number density 
in the energy interval {E,E + dE)), q is the source term, Dxx is the spatial diffusion coefficient along the regular 
magnetic field lines, Vc is the velocity of the Galactic wind, Dpp is the coefficient of the diffusion in momentum space, p 
is the momentum loss rate, and r/ and Tr are the time scales for fragmentation loss and radioactive decay, respectively. 

The transport equation is solved numerically and assuming a cylindrical symmetry, with halo boundaries at disc ra- 
dius i? = 20 kpc and half-thickness Zh as described below. We exploited a modified version of the GALPROP code [3^ . 
The main modifications consist in introducing by input the spatial and spectral profiles of the DM source (computed 
within the DarkSUSYpackage [13]), and in including the possibility of a spatially varying diffusion coefficient. 

In the following, we mainly consider one-zone models with isotropic diffusion, which can be regarded as the most 
extensively tested models of the recent past (see, e.g., Ref. ^] for a review). 

Our approach is to perform self-consistent tests in the local region and the parameters in Eq.[T]are chosen to strictly 
reproduce the local directly-observed spectra of nuclei and electrons. 

The goal of the paper is to study the possibility of disentangling the diffuse signals originated from two different 
sources, CRs and DM, having different spatial distributions. The CR injection source is confined to the Galactic 
plane, while the DM profile has a spherical shape. The region with intermediate and large z is thus the best target 
for the analysis. The propagation reshuffles the distribution of the two populations of electrons (and thus IC and 
bremsstrahlung signals), and the 7-ray signal associated to the decays of CR pions. The scaling of the signal along 
the z-direction is affected by almost any quantity entering in the transport equation, such as the description of the 
diffusion, the wind velocity, the magnetic field structure, and the interstellar radiation field (ISRF) distribution. 
Moreover, it is dramatically sensitive to the height of the propagation halo, namely, to the boundary condition along 
the z-axis. 

We are not interested in performing a full scan of the propagation parameters space and estimating the corresponding 
uncertainties in the CR spectra (see, e.g., Refs. ^E[13); rather, we want to investigate how different scalings along 
the z-direction due to different propagation models can affect the predictions for the signal to background ratio. We 
consider six benchmark scenarios of propagation and injection spectra, which are summarized in Table 1. In the 
following, we motivate our selection. 

Halo height: In addition to the "conventional" model having Zh =4 kpc (named BO), we consider two models 
of propagation in which the halo height has been set to Zh=^ kpc (model Bl) and z;i=10 kpc (model B2). The 
strongest constraints on the halo height is given by the "radioactive clocks", namely, unstable secondaries. Indeed, 
the ratio between stable and decaying isotopes is sensitive to the CR confinement time, which is in turn related to 
the halo height (and the diffusion coefficient). At present, the most precise measurements is the ratio ^°Be/^Be, with 
the unstable ^'^Be decaying in 10^ years. In Fig. [T^, we show the local interstellar spectra (LIS) and the modulated 
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spectra of the ^°Be/^Be ratio. The solar modulation is computed in the force field approximation [ssj. As expected, 
=4 kpc seems to be preferred by data. The model B2 is fully consistent with data at low energy (which are the 
most reliable), while it predicts a slightly smaller ratio than the two points at 1 GeV (which, however, have very 
large error bars). The model Bl shows, on the other hand, friction with data. Indeed, in the case with z^—l kpc, any 
model for diffusion and convection that can fit the ^"Be/^Be data dramatically increases the time spent by particles 
at low energy in the system, in order to let the ^°Be decay. This leads to an overproduction of secondary, at odds, 
e.g., with the B/C data at low energy. On the other hand, the properties of the local environment are highly relevant 
for radioactive species, since they typically travel a short distance before decaying. The presence of a Local Bubble 
surrounding the Sun can significantly increase the amount of decays (lowering the ratio) 39] . A precise estimate of 
this efi^ect, however, would require a detailed description of the local environment (probably, a collection of clouds) 
and this subject has not been fully addressed yet [3a |. For this reason, we consider the model with zi^=l kpc as 
disfavoured, but not definitively ruled out, and we do not disregard it. 

Convection: CR data at energies below the "knee" are in fair agreement with the description of the CR transport 
in the Galaxy through a difi^usive model, with just the relevance of convection in some regions under debate. Galactic 
winds are driven by active galactic nuclei or starbursts [40*1. Although the Milky Way have a star-formation rate 
significantly lower than starburst galaxies, the diffuse soft X-ray emission in the inner region of the Galaxy can be 
well explained by a kpc-scale wind [41]. This has motivated a renewed interest in model of cosmic-ray transport in 
which the convection plays a significant role |42|. The behaviour of the ratio of secondary to primary abundances 
disfavours a very high wind speed on the disc [33|. Moreover, in the outer region of the Milky- Way, there are no 
direct observational evidence of a gas-loss in a wind as for the inner region. On the other hand, it has been shown 
that models with a significant convection but an anisotropic diffusion or with two-zone of propagation [i^] (a 
purely diffusive zone close to the disc and the diffusive/convective zone outside), or with a wind speed increasing with 
distance from the disc [s^ , are consistent with CR data. 

In the CR models we are going to consider (i.e., one-zone of propagation and homogeneous diffusion), if the 
convection takes over at low energies, the secondary to primary ratio generally flattens, leading to conflict with data. 
The convection velocity Vc is thus very constrained on the plane, and we take Vc{z = 0) = 0. Models for thermally and 
cosmic-ray driven winds predict speeds increasing with the distance from the disc, with a linear scaling at intermediate 
z [i^. Radioactive isotopes and B/C data lead to the constraint: dvc/dz < 10 km s^^ kpc~^ [s^. We will analyze 
the effect of convection on the diffuse emission in the benchmark model B3, considering = 10 • {z/1 kpc) km s~^. 
In such model, electrons at high latitudes (and low energy) can stream away due to the wind transport, reducing the 
related IC and bremsstrahlung signals. 

Spatial diffusion and Reacceleration: The scattering of cosmic-ray particles on random hydromagnetic waves 
leads to diffusion and stochastic acceleration of the particles. The process of diffusion has a resonant character 
and it is mainly driven by the energy density associated to the random component of the magnetic field at the the 
resonant wave number of the scattering, fc^es = l/'^gi where Vg is the gyro-radius of the particle. In the quasi-linear 
approximation of turbulence, the perpendicular diffusion is subdominant, and the diffusion takes place mainly along 
the magnetic field lines, with the the form of the diffusion coefficient ^32*]: Dxx{,T,p) = l/3rgVp{SB{kres)/B)~'^, where 
Vp is the particle velocity and 6B is the amplitude of the random magnetic field. Combining data on fluctuations of 
the thermal electron density, interstellar cloud density, and interstellar magnetic field [i^ , the spectrum of the energy 
density of the interstellar turbulence in the nearby region (< 1 kpc) of the Galaxy is found to be well described by 
a power law, /c^^+". The data are consistent with a Kolmogorov-like spectrum (a — 1/3) over a large range (10® to 
10^° cm) of wavelength, and the inferred outer scale of turbulence is L~100 pc. Strong fluctuations on such large scale 
tend to isotropize the CR distribution. The diffusion coefficient can be then described by the form: D^x — PpDo Rqv 
with R being the rigidity of the particle, Pp = Vp/c and typical value for the diffusion coefficient in the Galaxy are 
Dq ~ 10^^ — 10"^°cm^s~^ [11]. Homogeneous diffusion can be regarded as a good approximation for the analysis of 
the nearby region. The picture at, e.g., large z is, on the other hand, more uncertain and spatially varying diffusion 
coefficient can quite significantly affect the predictions for the 7-ray signals in such region of the sky [43 , while local 
spectra were not highly modified, being spallation mainly confined in the Galactic disc. We consider such possibility, 
by describing the diffusion coefficient as D^x — Dxx exp(0/zs), with = 4 kpc, in the benchmark model B5. 

The stochastic reacceleration process refers to a diffusion in momentum which can be described by the coefficient 
Dpp — p^v^/{9 Dxx), where Va is the Alfven velocity associated to the propagation of the hydromagnetic waves [sl]. 
The expressions considered for Dxx and Dpp can describe the isotropic distribution at large scale of the Alfven and 
fast magnetosonic waves (interactions between particles and slow magnetosonic waves are usually small). Numerical 
and analytical treatments of the Alfven waves cascades typically leads to a Kolmogorov spectrum for the interstellar 
turbulence [1^; this case, due to the fast cascade rate, is not significantly affected by damping on cosmic rays and we 
disregard this effect. Typical values for the Alfven speed in the Galaxy are Ua ~ 10 — 30 km/s. 

In addition to the Kolmogorov case, we consider a Kraichnan spectrum of interstellar turbulence, i.e., a = 0.5. In 
this case, damping on cosmic rays can significantly affect the reacceleration mechanism and we include it following 
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the line of [43] (taking the damping constant g = 0.06 and the maximum free path length of 15 pc). This model of 
propagation is labeled B4 in Table 1. 

The diffusion coefficient Dq and Alfven velocity Va are tuned in order to reproduce the B/C ratio in all the 
benchmark models. Values are reported in Table 1. Spectra of the B/C ratio are shown in Fig. [TJj. Note that all the 
benchmark models are satisfactory in fitting the B/C data. 

In the rest of the paper, we are mostly sensitive to the high energy (> 10 GeV) description of the propagation. 
We check the reliability of our models, by performing a x^-analysis and comparing the predicted B/C ratio with data 
from the most accurate surveys, namely, CREAM ^491], ATIC [50,], HEA03 ^51], and CRN [5^], at > 3 GeV. Results 
are reported in the last column of Table 1 . 

Energy loss: All the energy losses are computed within the Galprop code as described in Ref. [ssj . 

Updated calculations of the ISRF [53j have estimated a quite different emission in the inner region of the Galaxy, 
depending on the assumptions on the metallicity gradient. The picture in the outer region is, however, basically 
unchanged, and our results can be only very mildly affected. 

The large-scale structure properties of the magnetic field are not extremely important as the turbulence properties 
in determining the diffusion. The strength is, on the other hand, crucial for the estimates of both synchrotron radiative 
emission and energy loss. It turns out that the latter is a subdominant component of the energy loss term of Eq. [T]in 
most regions of the Galaxy. A precise estimate of the magnetic field strength is thus marginally relevant for X- and 
7-ray emissions, while it becomes obviously very important when discussing radio and infrared signals. 

The magnetic field strength can be estimated from pulsar data as 5^]: B = Boexp{— ^j^^^" ), with Rg being the 
Sun-Galactic Center (GC) distance, Bq = 2.1 ± 0.3 and Rb = 8.5 ± 4.7 kpc (similar results from extragalactic 
sources ^55"]). Note that polarization observations refer only to the line-of-sight component of the magnetic field. 
Radio synchrotron measurements suggest higher values for the strength of the total field B, namely, 6 fiG near the 
Sun and about 10 /iG in the inner Galaxy (outside the GC), assuming equipartition between the energy densities of 
magnetic fields and cosmic rays |56] (this result is fairly in agreement with observations through the Zeeman splitting 
of atomic and molecular lines [57]). The radial scale length of the equipartition field is of about 12 kpc. On the 
other hand, analysis of the WMAP synchrotron foreground data (plus some assumptions on the CR distribution and 
turbulence model) can lead to [s^ Bq = 3/xG, Rb = 11 kpc, and Bturb/Bo = 0.57, not far from the estimate though 
rotation measures of pulsars. 

We consider the benchmark case B = 5 ■ exp[— (i? — i?s)/10kpc — ]z]/2kpc] fiG, with Rs = 8.5 kpc. 



6 



III. SOURCE TERMS 
A. Standard primary cosmic ray components 

There are strong indications that the main mechanism of acceleration for primary Galactic CRs, up to energies of 
100 TeV or so, is the scattering of CR particles with the strong shock wave fronts produced by supernova remnants 
(SNRs) in the circumstellar medium [59j. We will assume the primary CR source to be in the form: 



Qf (i?, z, E) cx i?"' cxp exp (^-^^ E'^- 



(2) 



where — 2.35, the radial length scale Rs — 1.528 kpc, and the vertical cutoff Zg = 0.2 kpc confines the source 
distribution to the Galactic plane. Neglecting discreteness and time variation effects, which could be eventually 
considered in connection to young nearby SNRs, the spatial part of the source function follows the mean SNR 
distribution in the Galaxy as derived from radio pulsar population surveys [63 |. This functional form, which is 
also in rough agreement with the gas distribution, traces the Galactic Type II supernova distribution and is highly 
suppressed in the Galactic bulge region. A contribution to primary cosmic rays from Galactic Type la supernova could 
be potentially relevant in the inner Galaxy (as one finds extrapolating from the distribution of old disc stars [61]), while 
it's certainly very subdominant in the local environment; since we will be mainly concerned with local observables, 
we are not sensitive to this this contribution and simply disregard it. 

Regarding the energy dependence in the primary CR source function, the theory of first order Fermi acceleration 
at astrophysical shocks predicts, for relativistic particles, a power-law spectral behavior, with the injection spectral 
index pmj ~ 2 in the limit of strong shocks [64I. Assuming that the local measurements are representative of the 
primary CR densities at the Sun galactocentric radius, we deduce the injection indices as well as source function 
normalizations matching the primary spectra after propagation to the local measurements. For cosmic-ray nuclei, at 
high energy, spatial diffusion is the dominant term in the propagation equation Eq. ([1]); the injection index can then 
be derived simply subtracting the scaling of the spatial diffusion coefRcient a from the spectral index of the locally 
measured flux, (3nuc — 2.7, namely f3inj,nuc — Pnuc — a, with a = 1/3 for a Komogorov diffusion, or a = 1/2 for a 
Kraichnan spectrum. On the other hand, energy losses leads typically to a softer spectrum for high energetic primary 
electrons 63]. The most precise measurement of the -I- e~ spectrum at high energy has been recently reported 
by the FERMI collaboration [ii]. The locally observed scahng with fie ^ 3.05, can be reproduced assuming 

Pinj,e — 2.35. At enrgy above 1 TeV, HESS data [H] suggest the presence of an exponential cutoff. We multiply 
Eq.[2]by exp(— (iJ/5TeV)^), which is found to reproduce well the experimental data. Finally, for both primary nuclei 
and electrons, we allow for a low energy break in the spectral indices, as suggested by the diffuse synchrotron and 
soft 7-ray emissions of the Galaxy [6^. 

The local fluxes, — Ni ■ Ef Vi/ {Att), for protons and primary and secondary electrons as computed for the six 
benchmark models selected in the previous Section are shown in Fig. O The corresponding injection spectral indices 
are reported in Table 1. 

Several authors have questioned the assumption that local measurements are representative of the CR density in 
the Galaxy, see e.g. [67| . Substantially different spectral indices at high energy or normalizations have been sometimes 
invoked for electron and/or proton spectra, e.g. to address [20] the excess above 1 GeV reported by EGRET in the 
observations of the diffuse gamma-ray emission from the Galaxy pJl] . This scenario seems now not strongly supported 
by preliminary data from the FERMI gamma-ray telescope. The spectrum of electrons with energy E >^ \ TeV is 



highly affected by the local picture, since they travel short distances due to the energy losses. Recent HESS data [65 1 
suggest the existence of at least one source of CR electrons within ^ 1 kpc. On the other hand, the steepening in 
the energy spectrum without a clear hardening before seems to exclude strong contributions from local sources. A 
different spectrum for protons is even more unlikely, since the scale over which protons diffuse is much larger than 
for high energy electrons, and their density is not sensibly affected by a spatially inhomogeneous energy loss term. 
We will not consider further this possibility; all results below are derived under the hypothesis that primary CR 
components can be normalized to the locally measured fluxes. 

B. Additional astrophysical sources for electron/positron cosmic ray components 

Secondary electrons and positrons derive from the decays of charged pions produced in the interaction of primary 
cosmic rays with the ISM along their propagation in the Galaxy. Their spectral index at the source is equal to the 
spectral index of primaries after propagation , i.e. close to Pnuc — 2.7, which is larger than the injection spectral 
index for electrons, Pmj^e — 2.35. The ratio of secondary to primary electrons is thus expected to decrease as the 
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FIG. 2: Local proton {Left Panel) and electron {Right Panel) spectra. Data are compared to the benchmark models of 
propagation BO (blue), Bl (green), B2 (red), B3 (orange), B4 (cyan), and B5 (magenta), described in Table 1. 



energy increases. To reverse this trend, and fit the sharp raise in the positron fraction detected by PAMELA |j| above 
10 GeV and up to 100 GeV, it seems unavoidable to introduce an extra electron/positron source with harder spectrum. 
Several possibilities have been discussed. Pulsars are well-motivated astrophysical candidates, see, e.g., 0, H, We 
consider here the alternative proposal of a possible production of secondary j e~ from hadronic interaction within 
CR sources. Being injected in the same region where cosmic rays are accelerated, they would have a harder spectrum 
than secondaries in the ISM, leading to a secondary to primary ratio increasing with energy jTo| . 

We model this additional component assuming that the spectrum at the sources is described by a power-law plus an 
exponential cutoff: E'^^^^-""-" • exp(— £'/£'c). We consider a Bohm diffusion inside the source, which implies a spectral 
index f3inj^sas — Pinj,nuc — 1 (^t high energy). The spatial part of this extra source is assumed to be the same as for 
standard primary components. The normalization follows instead from the requirement that PAMELA data can be 
fitted when including this additional term. 

Although the physical insight for this picture is different from the case in which the enhancement in the positron 
fraction is due to one or few nearby pulsars, from the point of view of testing the scenario through radiative emission 
the two cases are hardly distinguishable. In both cases local sources dominate the signal, and in both cases these 
sources are confined in the thin vertical layer where standard primary sources are confined. The discussion we present 
below for secondaries at the sources is then readily extendable to the pulsar scenario. 



C. Component from dark matter annihilations or decays 

A further possibility is that the extra component needed to explain the rise in the positron fraction is an exotic 
term due to dark matter in the Galactic halo. There are two possibilities: WIMP dark matter particles are stable 
but can annihilate in pairs injecting a given species i; the source term associated to this process is given by: 

Q1{r,E)^{aav)^x^{E), (3) 

where p{r) is the Milky Way halo mass density profile, assumed for simplicity to depend only on the spherical 
coordinate r, the mass of the dark matter particle, av the pair annihilation rate for typical velocities of dark 
matter particles in the Galactic halo (namely, in the zero temperature limit, as opposed to the finite temperature 
regime which applies in the early Universe), and dN^ /dE{E) is the number of particles i emitted per annihilation in 
the energy interval {E, E + dE). The second possibility is that dark matter particles have a long but finite lifetime, and 
the species i is injected in dark matter decays (for the interpretation of the PAMELA anomaly in terms of decaying 
DM, see, e.g., [sl,!^!, [t^jIzH); in this case the source function takes the form: 
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where Td is the decay rate and dNf/dE{E) is the number of particles i emitted per decay in {E^ E + dE). 

The distribution of dark matter in the Galaxy is rather poorly known, and one has to rely on large extrapolations. 
One possibility is to take N-body simulations of hierarchical clustering in cold dark matter cosmologies as a guideline. 
Numerical results indicate that dark matter halos can be described by density profiles that are sharply enhanced 
towards the Galactic center; there is still an on-going debate regarding how cuspy the profiles are, while, from the 
observational point of view, a tension has been often claimed between the profiles found in the numerical simulations 
and the dark matter density distribution as inferred from circular velocity or velocity dispersion maps, especially in 
small or low brightness objects. For what concerns our analysis, introducing a dark matter model with large density 
towards the center of the Milky Way would trigger discussion on signals generated in the central portion of the Galaxy. 
In line with other recent analysis, e.g., Refs. jl^, [13, [H, Ull, [T^I , strong constraints would arise. However these limits 
are strongly model dependents, relying on several extrapolations with respect to quantities that are measured and or 
tested with local observables, first of all the one on the dark matter density profile itself. We take here a different 
route and try to sketch results that depend as weakly as possible on extrapolations of local quantities. For the dark 
matter profile we assume a functional form with a large core radius in the central region of the Galaxy, the Burkert 
profile [73l |: 

^'^'") = 7n / ^ fi" I / / ^2^ ' 

(1 + r/rc) (1 + [r/rcY) 

where r is the distance from the GC; profile normalization po = 0.84 GeV cm~^ and the core radius = 11.7 kpc, 
corresponding to a local halo density of p(ro) = 0.34 GeV cm^^, are obtained for a model of the Milky Way halo 
with virial mass and concentration parameter being, respectively, Myij. = 1.3 x IO^^Mq and Cmr = 16. This model 
matches a whole set of available dynamical informations [t^], including, among others, constraints from the motion 
of stars in the Sun's neighbourhood, total mass estimates following from the motion of the outer satellites, and the 
Milky Way rotation curve. The Burkert profile, which was originally introduced as a phenomcnologically model from 
dynamical observations, from the theoretical point of view has been discussed as the equilibrium configurations from 
cold dark matter models in case the baryon infall happens with a sensible exchanged of angular momentum between 
baryons and dark matter [TSj . 

Note that for the dark matter density profile in Eq. ^ the gradient of the density profile in the region where most 
of the signal we will consider originates is very modest. Hence, in practice, the spatial signature of the annihilation 
source function, due to the scaling of the source function Eq. ([3|) with the number density of WIMP pairs, i.e. with 
/9^(r), is hardly distinguishable from that of the decay source function, scaling simply with the number density of 
dark matter particles, i.e. with p{r). 

Unresolved Galactic subhalos could, in principle, sligthly modify our conclusions, depending on their description. 
In particular, if the spatial distribution of subhalos is antibiased with respect to the host halo mass distribution 
(as found, e.g., in the Via Lactea simulation [76|). the emission at mid-high latitudes would be enhanced and the 
conclusions would be strengthened (while they would be weakened in the opposite picture). On the other hand, the 
estimates of mass function, spatial distribution and concentration of subhalos have great uncertainties and, taking a 
conservative approach, we disregard the contribution of substructures to the total DM density. 

If the dark matter particle mass is not known, it will be difficult to discriminate the two cases even from the 
spectral shape of the sources. To reproduce the PAMELA positron rise, a positron dark matter source with hard 
spectrum is required; this is possible in case of prompt annihilation or decays into leptons or weak gauge bosons. The 
second possibility is strongly constrained by available measurements of the antiproton cosmic ray flux (see, e.g., [l^). 
since and are copious sources of antiprotons. We will consider annihilations in pairs into monochromatic 
e~^/e~ {Ef. ~ M^) and monochromatic e+/e~ from a two-body final state in the decay {E^ — M^/2), as well 
cascades initiated by annihilations or decays into and /t^ . In the cases with e^/e" and as final 

states of annihilation/decay, the resulting /e~ source comes together with a gamma-ray component stemming from 
radiative processes at emission; we consider the model independent part of the final state radiation (FSR) terms (limit 
of ^ me) and neglect eventual model depend terms for specific WIMP pair annihilation cases (sometimes referred 
to as virtual internal bremsstrahlun g |77l . [78j ) or for the case of a DM candidate annihilating into new light particles 
which in turn generate e~'"/e~, e.g., [79| (in this last case the FSR component may be reduced with respect to what 
we are implementing in the rest of the paper). For the case mediated by r+/T~, on top of the FRS term, there is 
a substantial 7-ray flux induced by neutral pion production and their subsequent decay into 2 photons; this term is 
accounted for linking to the Pythia Montecarlo simulations as provided by the DarkSUSY package. 

The 7-ray signal is directly related to the injection source of Eqs. [3] and 21 The computation of the radiative 
emissions requires the solution of Eq. [Ijwith the e~^/e~ injection source described in Eqs. [3] and [4] We perform the 
calculation of this term with the Galprop package, analogously to the previous contributions. Relevant formulas to 
infer some scalings of the results (in case propagation is treated in an approximate form) are described, e.g., in Sec. 3 
of Ref. [si. 



9 



FERMI 
ATIC-1,2 
PPB-BETS 
BETS 
HEAT 
HESS 

HESS syst. band 




10' 10" 
E[MeVl 



0.01 



PAMELA data 

secondary at source 

DMe 

DM\l 

DM-u 




secondary in ISM 



A... 

; \ 
•. \ 
■■ \ 



10 



10" 
E [GeVl 



10" 



FIG. 3: Left Panel: Local electron+positron spectrum in the "conventional" model BO. We show primary CR electrons (blue 
solid), secondary CR electrons + positrons produced at the source (black dotted), and electrons + positrons induced by WIMP 
annihilation in the model DMe (blue dotted), DMr (blue dashed), and DM^ (blue dashed-dotted). Right Panel: Local positron 
fraction above 10 GeV in the "conventional" model BO. Best-fit curves including the contribution of the positron sources shown 
in the left panel are plotted with thick lines. They are obtained by allowing mild variations in the normalization and spectral 
index of primary and secondary electrons/positrons spectra, as described in the text. Contribution from secondary positrons 
produced in the ISM is shown with thin lines. Line-styles as in the left panel. 



IV. SIGNAL VS BACKGROUND 



A. Electron/Positron spectrum 



For sake of clearness, we concentrate our analysis on three benchmark WIMP scenarios, named DMe, DMr, and 
DM^, and summarized in Table 2. Injecting monochromatic e'^ /e~ yields, the model DMe induces a very steep 
electron/positron spectrum. This is true also for the model DM/x, while for the model DMr, which assumes t'^ /t~ 
are produced, the spectrum is quite hard but significantly smoother. A dark matter scenario fitting the PAMELA 
positron excess typically lies in between these three cases. Note that in the case of DM annihilations/decays producing 
such final states after a number of steps mediated by new exotic light particles, the final spectrum of j e~ is softened 
with respect to the picture considered in the following (to an extent depending on the number of steps and on the 
properties of the new particles). 

We start with the case of pair annihilating WIMPs. The explanation of the PAMELA data requires a WIMP 
mass > 80 GeV. WIMPs with < 1 TeV and inducing a very hard spectrum of /e~ are disfavoured by 
the FERMI data on e+ -I- e^ (assuming a power-law for the injection spectra of primary electrons). We consider the 
scenario of a WIMP annihilating directly into e^/e" as a toy benchmark case to show that for any DM model inducing 
a hard spectra of e^/e" fitting the PAMELA positron fraction, the associated diffuse emission is at a detectable level, 
independently on the mass (for > 150 GeV). A heavy WIMP (i.e., > few hundreds of GeV) with the same 
final state would actually induce an IC emission which is already strongly constrained by EGRET data [3l| . For the 
DMe benchmark scenario, we consider a WIMP mass of 300 GeV. A fight (i.e., < 300 GeV) WIMP candidate 
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TABLE 11: WIMP benchmark models. The annihilation rates are listed in the case of the "conventional" propagation model 
BO. 
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FIG. 4: Electron vertical profile at i? = 8 kpc and E — 200 GeV. Left Panel: Propagation model BO: We show primary CR 
electrons (blue solid), secondary CR + e~ produced in the ISM (blue short dashed), secondary CR e"^ + e~ produced at the 
source (black dotted), and + e~ induced by DM annihilation in the model DMe (blue thick dotted). For comparison, we 
plot the distribution of a source scaling as pjjM (black solid) , with an arbitrary normalization. Central Panel: The same of the 
left panel, but adding the propagation models Bl (green) and B2 (red) and considering only the contributions from primary 
electrons and DM induced electrons + positrons. Right Panel: The same of the central panel, but in the propagation models 
B3 (orange), B4 (cyan), and B5 (magenta). 



annihilating into r+r" does not provide an explanation to the PAMELA data. The 7-ray emission induced by a 
heavy (i.e., > 500 GeV) DM candidate with the same annihilation modes (and fitting the PAMELA excess) could 
violate constraints by ACT measurements of dwarf satellites [8lj and of the Galactic Ridge [t^]. We refer to a sample 
case standing in middle of this mass interval, namely a WIMP mass of 400 GeV (DMr). A WIMP annihilating into 
with mass > 1 TeV can give a very good fit to, simultaneously, the PAMELA and FERMI data. In the 
benchmark case DM/i, we consider M-^ = 1.5 TeV. 

We tune the annihilation cross section in order to fit the positron fraction measured by PAMELA above 10 GeV 
and the e"*" + e~ local spectrum measured by FERMI. The fit is performed adding the dark matter component to 
the background contribution from standard primary CRs and the secondaries produced along propagation, computed 
in the same propagation model configuration. We allow the normalization and spectral index of both primary and 
secondary CR electrons to vary with respect to the cases reported in Table 1 (within the ranges A(3inj,e = ±0.5 and 
^f3inj,nuc = ±0.1), reflecting some uncertainties on their modeUing. Within these assumptions, the models describing 
the DM component with the benchmark scenario DMe cannot reproduce simultaneously the two sets of data. They 
would require some exotic assumptions for the primary injection spectra. Therefore, in this case we restrict the analysis 
to the fit of the PAMELA dataset only, requiring the e+ + e~ spectrum induced by DMe alone to be consistent with 
the FERMI observations. 

The best fit values for the DMe scenario are: CTqU = (2.5,2.3,2.7,2.5,2.7,2.1) • 10"2'*cm3s"\ for DMr scenario 
are: aaV = (6.6, 4.2, 6.1, 8.0, 6.6, 5.4) • 10-^4 cmSg-i, ^^^^ j^y^^ scenario are: aaV = (2.5,2.0,2.3,2.8,2.8,2.4)- 
10~^^cm^s^^ in the (BO, Bl, B2, B3, B4, B5) propagation model, respectively. For the DMr and DM/i scenarios, all 
the models provide a good fit to the data. The theoretical curves do not reproduce the data in an extremely precise 
way only for what concerns the positron fraction at energies below 20 GeV in the propagation models Bl and B3. 
In such cases, the production of secondary positrons from primary CR is enhanced (see Fig. ^jp) and the background 
alone tends to be in conflict with PAMELA data at low energies. We plot in Fig. [3] the obtained e+ + e^ spectrum 
and positron fraction of the three benchmark DM scenarios in the propagation model BO. 

Note that, as expected, the best-fit values for the annihilation rates are significantly larger than the thermally 
averaged annihilation rate at the time of decoupling of WIMPs, namely, aaV = 3 • 10~^^ cm^s~^. Many explanations 
have been proposed to motivate such mismatch, like, e.g., non thermal production mechanism, non-standard cosmology 
at the time of DM decoupling, and Sommerfeld enhancement. 

As we already mentioned, in the case of Burkert DM profile, the differences in the spatial distribution of a signal 
scaling with p^, as in the WIMP case, or with p, as in the decaying DM case, are very mild. Therefore, our results can 
be rephrased in term of benchmark decaying DM scenarios with decaying modes analogous to the annihilation modes 
introduced above, and with mass scale doubled. The corresponding decay rate can be approximately estimated by 
T ~ (CTaw)"^ • M^/po ~ 3.3 • 1026(^,^^/;L00GeV) • (10-^4 cm3s-V(CTat;)), where is the mass of the WIMP in the 
annihilating DM scenario. 

As discussed in the previous Section, for comparison, we consider also the case of an extra non-standard component 
due to the production of secondary e+/e^ inside CR sources. In Fig. [3l we plot the e+ -I- e^ spectrum and positron 




FIG. 5: 7-ray diffuse spectrum at intermediate latitudes (fO° < b < 20°), integrated over longitudes 0° < Z < 360° and compared 
to the FERMI preliminary data [2f|]. Left Panel: Emission in the propagation model BO. The CR (primary+secondary) spectra 
associated to yr^-decay, IC, and bremsstrahlung are shown by thin dotted lines. The thick solid blue line is the sum of the three 
components. The solid black line shows the extragalactic background in the model described in the text (thick) and fitted from 
EGRET data ^ (thin). The IC + FSR emission associated to the WIMPs DMe and DM/x are shown by thick dotted and 
thick dashed-dotted lines, respectively. The IC + 7-ray from 7r°-decay signals induced by the WIMP DMr are shown by thick 
dashed lines. Central Panel: Emission in the propagation models Bl (green) and B2 (red). Same line styles of the left panel. 
Right Panel: The same of central panel, but for the propagation models B3 (orange), B4 (cyan), and B5 (magenta). 



fraction in the propagation model BO for a population injected with Pinj,sas = 1-36, an energy-cutoff ai Ec = 1 TeV 
and normahzation tuned to the best fit value for the PAMELA and FERMI data (again allowing mild variations in 
the normalization and spectral index of primary and secondary electrons) . 

The absence of sharp features in the e+ -|- e~ spectrum detected by FERMI implies that the contribution of the 
source of positrons responsible for the PAMELA excess to the total e+ + local flux or to the associated diffuse 
emission, can be hardly singled out. If the spatial distribution of the sources of such "exotic" components traces the 
CR sources, the last sentence would be true everywhere in the Galaxy. 

In Fig. m we plot the vertical profiles of the electron number density distributions, at the local radial distance 
R = 8 kpc and E = 200 GeV. For this slice of the Galaxy, at rather large distance from the GC, the determination 
of propagation model parameters as derived by matching the LIS of nuclei is rather robust. We are focusing on 
some typical energy at which electron and positron sources relevant for the raise in the positron fraction are also a 
significant contribution to the total population of -\-e~ . In Fig.|4^ we consider the propagation model BO, and plot 
the vertical profile of CR primary electrons (solid), secondary e"*" + e~ produced in the ISM (short-dashed), secondary 

+ e~ injected at the source (dashed-dotted), and e"*" -|-e~ flux induced by WIMP annihilations in the DMe scenario 
(dotted). All these cases but the latter follow a distribution which is mainly confined to the disc (although broadened 
by the diffusion). The DM- induced component is instead much flatter (we plot for comparison the proflle of the DM 
injection source oc P^jm)- is the dominant component at intermediate and large z. We thus expect the associated 
radiative emission to dominate at intermediate and high latitudes. 

In order to understand how this conclusion is dependent on the propagation model considered, we show the cases 
of the propagation model Bl and B2 (plus again BO, for comparison) in Fig.[4]3, and of B3, B4, and B5 in Figllj:. In 
these figures, we do not plot the vertical profiles of secondary e'^ -t- e~. Their shapes are analogous to the CR primary 
electrons profile and the rescaling factor is roughly the same as in Fig. 3^. Note from Fig. [JJa that, as expected, as 
the boundary of propagation Zh increases (decreases), the region at which the DM-induced component is dominant 
becomes larger (smaller). From Fig.|3h, we conclude that at high energies the effect of convection (model B3) on the 
shape of the vertical profile is negligible. The same conclusion applies also to the effect of varying the spectral index 
of the diffusion coefficient (model B4). In the model B5, the population of electrons induced by DM at high z is 
mildly reduced with respect to the model B2 since the spatial diffusion coefficient is increasing with z, and electrons 
and positrons are less efficiently confined. 



B. 7-ray emission 



The discussion in the previous Section pointed out that, in order to detect a DM-induced signal in the diffuse 
emission of the Galaxy, intermediate and high latitudes are the best targets. 

At high latitudes, the diffuse extragalactic gamma-ray background (EGB) is expected to become the dominant 
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background component. To estimate the level of the extragalactic emission in the FERMI preliminary data (2l| 
reported in Fig. O we rely on the EGRET data and we consider the fit obtained in Rcf. 20] (upper black curve). 
The sharp increase in sensitivity to point sources of the FERMI telescope with respect to EGRET, may, on the other 
hand, lower significantly such term. In three months of observations, FERMI has already detected an amount of 
individually resolved active galactic nuclei (which are believed to be the main component of the EGB) corresponding 
to ~ 7% of the EGRET extragalactic diffuse gamma-ray background [m . We consider a model for the contribution 
of unresolved blazars as in Ref. [83] (lower black curve in Fig. O and [7]), estimating the FERMI point source sensitivity 
as 1.6 • 10~^cm~^s~^, roughly corresponding to 3 years of observations. Another crucial ingredient to estimate the 
diffuse extragalactic radiation is absorption of gamma-rays at high energies, mainly due to pair production on the 
extragalactic background light emitted by galaxies in the ultraviolet, optical and infrared bands. We consider the 
parametrization of this effect in Ref. [84] , as derived in the context of the ACDM cosmological model. 

In Fig. we plot the 7-ray diffuse spectrum at 10° < b < 20°, integrated over longitude (0° < ^ < 360°), and 
compared to the FERMI preliminary data. These measurement do not confirm the EGRET excess in the GeV 
energy range, with the level of the detected diffuse flux being significantly reduced. In Fig. [5^ we show the case 
of the "conventional" propagation model BO. In the previous Section, in order to determine the best fit values 
of the annihilation cross sections, we allowed the normalization and spectral index for primary and secondary CR 
electrons/positrons to vary with respect to the benchmark cases reported in Table 1 (see Fig. ^^). On the other 
hand, such mild variations in the spectra induce negligible variations in the associated diffuse emissions. Thus, in the 
following, for sake of clearness, we consider the background CR emission as given by the benchmark models of Table 
1. 

The first remark is that the sum (blue solid line) of three CR components (blue thin dotted lines), namely, IC, 
bremsstrahlung, and 7r°-decays, plus the extragalactic background contribution (black solid line), can approximately 
account for the measured flux at i? < 10 GeV (note that propagation models have not been tuned to do so, while we 
are just extrapolating from the LIS of nuclei). Exotic components, claimed in order to explain the EGRET excess, 
are now signiflcantly constrained, at least at mid-latitudes^. 

In the same plot one can see that the 7-ray flux induced by our benchmark DM models is more than one order 
of magnitude smaller than the detected flux at i? < 10 GeV, while it becomes comparable to or higher than the 
background a,t E > 100 GeV. At such energies, both the IC and FSR signals are relevant in the model DMe (thick 
dotted line) and DM/x (thick dashed-dotted line), while in the model DMr (thick dashed line) the flux is driven by 
the TT^-decay emission. Note that the in the DMe scenario the final e^/e~ spectrum is very sharp and thus the peaks 
of the IC emission on CMB and on infrared and optical starlight are clearly visible, while in the DM/x and DMr cases 
they are smoothed. 

The 7-ray diffuse spectrum at intermediate latitudes for the propagation models Bl and B2 is shown in Fig. [5Jd. 
As expected, the IC signal associated to the WIMP DMe and DMfi is enhanced (reduced) than the CR signal for 
the model B2 (Bl) with respect to the "conventional" model, making the WIMP scenario more easily detectable 
(undetectable). Moreover the CR secondary positron spectrum in the model B2 (Bl) is reduced (enhanced), see 
Fig. In the fit of the PAMELA excess, this leads to an enhancement (decrease) in the normalization of the DM- 
induced positron spectrum, and thus in the annihilation cross section, implying a higher (smaller) signal associated 
to TT^-decay and FSR. 

From Fig. [5}:, one can see that the level of CR diffuse spectra at high energy in the model B3, B4, and B5 is almost 
identical to the "conventional" case. The DM- induced emissions are also analogous. 

All the CR benchmark scenarios considered in the discussion induce emissions which are fair in agreement with 
the FERMI preliminary data. We find a slight mismatch only in the model Bl at high energy, caused by the lower 
level of IC and bremsstrahlung radiations as follows from the fact that we are introduced a smaller region for j e~ 
propagation and thus a smaller region for radiative emission, and in the model B4 at low energy, due to the choice of 
spectral indices (see the e^/e" spectrum at ~ GeV energies in Fig. [2b). 

As already emphasized a few times in the paper, the CR emission is almost confined on the disc (i.e., it is elongated 
in the i? direction and rapidly decreasing along z). Therefore, the diffuse emission in (6, Z) coordinates (namely, the 
fiux integrated along the line of sight labeled by the angles (&, /)) tends to have a fiatter longitudinal profile as the 
latitude increases. This is true also for the DM components (being the DM profile and the ISRF distribution, relevant 
for IC emission, decreasing with z), but to a smaller extent, especially at intermediate latitudes. We expect therefore 
larger ratios for the DM signal to CR background at longitudes close to the center of the Galaxy; on the other hand. 



^ Other observations reported by the FERMI LAT telescope (e.g., Vela pulsar 'ss'l) go in the same direction, namely, reporting a reduced 
flux at GeV energies with respect to the EGRET observations. The current most likely interpretation of the EGRET excess is thus an 
instrumental bias. This would imply that a significant contribution from exotic components at few GeV is severely constrained in any 
portion of the sky. 
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FIG. 6: Longitudinal profile of the 7-ray diffuse emission at intermediate latitudes (10° < fe < 20°) and E = 150 GeV. Left 
Panel: Emission in the propagation model BO. The CR profiles associated to 7r°-decay and IC are shown by thin solid lines. 
The solid black line shows the extragalactic background in the model described in the text. The profile of the IC emission 
associated to secondary e"*" + e~ produced at the source is shown by the black dotted line. The IC + FSR emissions associated 
to the WIMPs DMe and DM/x are shown by thick dotted and dashed-dotted lines, respectively. The IC + 7-ray from vr^-decay 
signals induced by the WIMP DMr are shown by thick dashed lines. Central Panel: Emission in the propagation models Bl 
(green) and B2 (red). Line styles as in the left panel. Right Panel: The same of central panel, but for the propagation models 
B3 (orange), B4 (cyan), and B5 (magenta). 

our procedure for estimating parameters in the propagation model is driven by local data, and thus our extrapolations 
are more reliable outside the inner region. 

In order to show that the detectability of a DM component inferred from Fig. [5] is not only related to the picture 
in the central region of the Galaxy, in Fig. [S] we plot the longitudinal profile for the diffuse emission at E = 150 
GeV and for 10° < b < 20° (namely, the same latitudes of Fig. [S]). Note that the diffuse component induced by DM 
dominates over a large range of longitudes for any benchmark scenarios. The picture is completely different with 
respect to the signal induced by secondary e"*" + e~ produced at the source, which is subdominant everywhere, having 
a spatial profile of injection identical to the primary sources. The CR emission associated to 7i"°-decay is shown only 
in the "conventional" case (Fig. [5^). At high energy, it is analogous for all the models of propagation, being mostly 
generated close to the disc. It is highly irregular due to the variation of the gas density. The longitudinal shape 
of the 7r°-decay signal in the southern hemisphere (—20° < b < —10°) would be different but with an overall size 
roughly analogous. From the figures, one can see that the longitudinal shapes of the signals are analogous for all the 
benchmark propagation models. 

In Figs. [7]and[8l we repeat the same analysis at higher latitudes, namely, 50° < b < 60°. As expected from the 
analysis of the previous Section, in this case we find a larger ratio between the signal induced by DM to the emission 
associated to CRs. On the other hand, the EGB starts to play a major role. The detection prospects for the WIMP 
scenario DMr are again very favourable in all the propagation models. The emission induced by DMe and DM/i is 
also detectable, being, roughly, of the same level of the sum of the backgrounds a.t E > 100 GeV. This is no longer true 
at higher latitudes, where the EGB takes over and such emission becomes too faint to give a clear signature. Fig. [SI 
shows that, as explained in the discussion above, the longitudinal profiles become flatter than at lower latitudes. The 
emissions come mostly from the local region and therefore these predictions can be assumed as rather robust. 

Note that the enhancement in the DM-induced IC emission in the propagation models with Zh = 10 kpc (B2 and 
B5) with respect to the "conventional" case (zh — 4 kpc) is more significant than at intermediate latitudes, and 
viceversa for the model Bl. The B2 case is more favourable than the B5 model; in the latter the e+/e~ population 
is slightly depleted at large z since the spatial diffusion coefficient increases in such region. The predictions in the 
models B3 and B4 are again analogous to the "conventional" case. 

The level we predict for diffuse 7-ray fluxes is about E^J ~ 1 — 3 ■ 10^"^ MeV cm~-^ s~-^ sr^^ at £' > 100 GeV (see 
Figs. OS]). Considering the FERMI performances stated in Ref. (roughly, an effective area of Aeff = 8 • 10^ cm^ 
and a field of view FoV = 2.4 sr), the expected number of counts, for an energy bin size of AE^ = 50 GeV, is about 
Nj > 70 sr^^ yr^^ . We deduce that the diffuse 7-ray spectra as predicted in Figs. [5] and [7] can be detected with 
a statistical error smaller than 10% in 1 year of observation. The precise description of longitudinal and latitudinal 
profiles requires, on the other hand, some years of observations. Combining different slices of the sky, however, the 
disentanglement between the CR source having a " disc" shape and the WIMP induced source having a spherical shape 
will be feasible in the forthcoming future. Full sky-maps, at 150 GeV for the Tr^-decay signal associated to primary 
CR and DMr, and for the IC emission associated to CR electrons and DM/x is shown in Figs. and [TUl Differences 
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FIG. 7: 7-ray diffuse spectrum at high latitudes (50° < b < 60°) integrated over longitudes 0° < I < 360°. Line styles and 
colors as in Fig. [5] 

in morphologies for the various components are indeed very clear. 

C. Radio and infrared emission 

Now we turn the discussion on the synchrotron emission in the radio and infrared bands. Electrons and positrons 
injected by DM or CR source interact with the Galactic magnetic field, giving raise to a synchrotron radiation. Due 
to the spectral behaviour, the synchrotron emission is the dominant component of the Galactic diffuse emission at 
low frequency. The sky-map of Ref. [U at 408 MHz is the standard calibration for the synchrotron diffuse signal 
(although it could include a significant amount of unresolved sources) . Foreground estimations in the WMAP data [88| 
suggest a spectral index for the synchrotron emission ~ 3, at frequency up to 60 GHz. (An anomalous component 
has be claimed to be present in the innermost region of the Galaxy, a result which depends on the template used for 
the foreground estimation. The associated spectral index turns out to be harder than 3. Such component, dubbed 
"WMAP haze", has been associated to be a possible DM signal due to WIMP annihilations [H, [iH, [O, [ill . Since 
the haze is associated to the central portion of the Galaxy, we will not discuss it here.) 

In Fig. [Til we show the emission associated to primary+secondary CR electrons in the " conventional" model at 
intermediate latitudes. Matching the diffuse emission induced by CRs with the observed synchrotron emission in the 
whole Galaxy is beyond the goal of this paper. Note, however, that the spectral index is very close to 3, as required. 
The overall normalization is also very close to the one estimated by the WMAP team. 

Again, in order to explore a possible DM signal, the region at intermediate and large latitudes is the best tar- 
get. Indeed, the magnetic field slowly decreases outside the disc (we adopt the benchmark case B = 5 exp[— (i? — 
i?o)/10kpc — |2;|/2 kpc] /xG, as described in Section [Hi) , and, thus, the signal induced by a spherical DM profile is 
non-negligible at intermediate and large z. 

The peak of the synchrotron emissivity is approximately at the frequency i^p ~ 4.7MHz • {Ee/1 GeV)^ /iG), 
with Ef, the electron energy and B the ambient magnetic field. The injection of high energetic electrons {E^, > 
100 GeV) induce a signal peaked at frequencies > 100 GHz. It follows that this represents the best frequency range 
for the search of a synchrotron emission induced by a DM candidate injecting an electron/positron yield at the level 
of the PAMELA excess. Such frequencies are above the WMAP range, while it could be possible to search for such 
a component with forthcoming PLANCK satellite, which will have detection channels up to a frequency of 850 GHz. 
We show in Fig. [Til the signals induced by the benchmark DM candidates introduced above, as computed in the 
propagation model BO. The detection prospects are less favourable than in the 7-ray band. Only the benchmark 
scenario DMe can induce a significant component at intermediate latitudes. 

Moreover, an extra uncertainty in the analysis at infrared frequencies is given by the fact that thermal emission 
rather than synchrotron is expected to dominate the foreground. Focusing the analysis on the spatial distribution 
and on polarization data could, however, help to disentangle a DM-induced synchrotron signal. A firmer statement 
on such a possibility with the PLANCK satellite deserves a more detailed dedicated study. 

In principle, data in the X-ray frequency band can also be used to test the diffuse continuum emission induced by 
DM pair annihilations or decays. On the other hand, the fluxes detected by INTEGRAL SPI [89^ and COMPTEL [9q| 
in the inner region of the Galaxy is a factor of few above the estimated CR diffuse emission (as computed in the 
"conventional" model). It follows that, generically, DM-induced fluxes are in turn less constrained by these data-sets 




FIG. 8: Longitudinal profile of the 7-ray diffuse emission at higli latitudes (50° < b < 60°) and E — 150 GeV. Line styles and 
colors as in Fig. [S] 



with respect to the 7-ray and radio observations. Moreover, the case of a DM candidates which offers an interpretation 
of the PAMELA excess is even harder to be detected, since the hard injection spectrum of the DM-induced e^/e" 
leads to a signal associated to the IC scattering with the CMB component which is subdominant with respect to the 
CR contribution, having the latter a softer spectrum (see Figs [5] and [T]) . 



V. DISCUSSION AND CONCLUSIONS 



We have presented a study of the diffuse 7-ray emission induced by particle-DM models which can account for 
the PAMELA positron excess. The analysis has been focused on signals induced by pair annihilations (or decays) 
associated to the smooth component of the Milky Way DM halo, in a configuration with no significant enhancement 
in the DM density in the Galactic center region, hence with none of the DM signals usually discussed from the center 
of the Galaxy at a detectable level. The goal was to present an analysis based on minimal sets of assumptions and 
extrapolations with respect to locally testable or measurable quantities (for comparison, e.g., Ref. [T^I considers also 
DM profiles with a central cusp, finding strong, though heavily model-dependent, constraints from GC observations). 
In this respect, e.g., we are not directly comparing to the possibility of detecting DM fiuxes from Milky Way satellites 
or other dark matter dominated objects, avoiding the extrapolations that are required when addressing this kind of 
predictions. Estimates of extragalactic 7-ray background from unresolved DM structures lead typically to a flux below 
the astrophysical EGB (unless a substantial enhancement stems from populations of dense substructures), which is 
mainly related to unresolved blazars; this component is therefore marginally relevant for our analysis (although it could 
be possibly detected studying its anisotropy pattern [9l|). On the other hand, a contribution of unresolved Galactic 
subhalos to the diffuse emission of the Galaxy could be substantial Our conclusions would be strengthened, in 
particular, in case of a spatial distribution of subhalos antibiased with respect to the host halo mass distribution [76|, 
which would give a larger diffuse emission at high latitudes and longitudes. While not changing the general picture, 
the estimation of this component depends on several quantities which are not well-known, such as the mass function, 
the spatial distribution and the concentration of subhalos. 

The presence of a Galactic dark disc can in principle affect our conclusions if the DM density associated to the disc 
is much higher than the density of the halo. However, this is quite unlikely, as also found in the simulation of Ref. [o^ . 
which proposed the existence of a dark disc in Milky Way-sized galaxies. Moreover, we consider, for simplicity, a 
spherical shape for the DM halo, which is actually an approximation. Significant deviations from sphericity could 
induce a different scaling of the signal in the longitudinal direction with respect to the latitudinal direction. However, 
our general conclusions will be unchanged, being crucially based only on the assumption that the DM halo is extended 
along the z-direction well further the stellar disc. 

We have compared the DM signal to standard astrophysical components, analyzing spectral and spatial features to 
derive which are the best energy band and region of the sky to disentangle the DM contribution from the Galactic 
plus extragalactic background with the FERMI LAT telescope. Being the sources of acceleration for CR (namely, 
Supernovae Type II) distributed within a thin disc, the ratio of the associated background versus the DM signal 
decreases as the latitude increases. The extragalactic background has, on the other hand, a flat spatial distribution, 
and dominates at very large latitudes. Intermediate and moderately large latitudes (10° < b < 60°) are thus the most 
promising portion of the sky for the search of a diffuse emission induced by Galactic DM. 
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FIG. 9: Sky-map at 150 GeV of the emissions associated to Galactic primary+secondary CRs in the "conventional" model BO. 
The intensity is shown in logarithmic scale and units [MeV cm~^ s~^ sr~^]. Left Panel: Inverse Compton radiation. Right 
Panel: 7r°-decay emission. 
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FIG. 10: Sky-map at 150 GeV of the emissions induced by WIMP annihilations in the propagation model BO. The intensity is 
shown in logarithmic scale and units [MeV cm~^ s~^ sr~^]. Left Panel: Inverse Compton radiation in the DM/i scenario. Right 
Panel: Tr^-decay emission in the DMr scenario. 

Fits to the PAMELA positron excess require sources with a rather hard spectrum. The preferred DM-related 
interpretation is a scenario with prompt emission of leptonic final states; we have discussed in detail three sample 
benchmark cases, namely emission of monochromatic e~^/e~, the case of /t~ yields, giving a shallower e'^ /e~ 
spectrum, and the channel as final state of annihilation/decay. We have disregarded WIMP models with 

annihilation channels producing soft spectra, like quark-antiquark pairs, and with weak gauge boson final states, 
which would overproduce cosmic-ray antiprotons 94]. 

We find that for all the benchmark models, after tuning the source rate to match the level of the PAMELA positron 
excess, the DM-induced 7-ray flux becomes comparable to the background at energies E > 100 GeV. In case of pairs 
annihilation into monochromatic e^/e~ or into fi'^ / ^~ , 7-rays arises from IC scattering of the propagating e^/e~ 
and FSR processes at emission; the scenario is detectable for sufficiently heavy DM candidates (M^ > 150 GeV; 
the PAMELA positron data require > 80 GeV). In the case of a DM candidate annihilating into r+r~ with a 
mass of few hundreds of GeV, as required by the PAMELA excess without violating ACT constraints, the detectable 
7-ray component is due to the emission from Tr^-decays, which peaks at, roughly, one-third of the DM mass. These 
conclusions can be straightforwardly extended to a decaying DM scenario: the mass scale just needs to be doubled 
and the decay rate properly adjusted. In case of a DM candidate annihilating into a new light particle which in turn 
decays into leptons, the FSR is generally reduced (depending on the model). The total emission can be thus mildly 
reduced with respect to an analogous WIMP case, but remains still sizable. 

The statement regarding the detectability of the induced Galactic diffuse 7-ray flux has a marginal dependence 
on the model implemented to describe the propagation of charged particles in the Galaxy. We introduce a locally 
self-consistent picture, implementing for e^/e~ (as well as for primary CR protons and all other relevant components) 
a transport equation with propagation parameters tuned to reproduce the LIS of cosmic rays. In particular, being 
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FIG. 11: Synchrotron radio and infrared emissions for CR primary + secondary electrons in the "conventional" model BO (solid) 
and for the benchmark WIMP scenarios DMe (dotted), DMr (dashed), and DMfi (dashed-dotted) . Left Panel: Spectrum at 
mid-high latitudes (10° < b < 60°) integrated over longitudes 0° < Z < 360°. Right Panel: Latitudinal profile at "small" 
longitude {I — 10°) and at frequency f — 650 GHz. 

mostly sensitive to the picture at GeV to TeV energies, we fit propagation parameters to B/C data at i? > 3 GeV. 
We consider as a reference model the " conventional" model introduced by the GALPROP collaboration [s^ . We then 
discuss the variation of the vertical profiles of CR and DM-induced electrons and of the associated diffuse emissions 
by changing the halo boundaries of propagation, introducing an advection term, modifying the spectral index for 
diffusion, and considering a spatially varying diffusion coefficient. We have found that the signal to background ratio 
is enhanced if the halo boundaries in the z-direction are extended with respect to the "conventional" case {zh = 4 
kpc), and reduced if they are restricted. The latter possibility is, however, disfavoured by CR data. For all the other 
propagation models, the detection prospects are analogous to the "conventional" case. 

We have also discussed the possibility of detecting radiative emissions at other wavelengths. In the radio and infrared 
bands, the synchrotron component associated to a scenario in which DM injects a very hard e"*" + e~ spectrum (like, 
e.g., a WIMP with a large branching ratio of annihilation into e^/e~) is expected to be a significant component of 
the synchrotron diffuse emission of the Galaxy at few hundreds of GHz. This frequency range can be probed by 
the forthcoming PLANCK mission. On the other hand, this picture is strongly constrained by the FERMI data on 
e"*" + e~, and the detection prospects at radio frequencies are, typically, less favourable than in the 7-ray band. For 
the class of DM candidates explaining the raise in the positron fraction, the X-ray energy band is even less favourable, 
being the IC spectrum induced by DM harder than the one associated to CRs. 

To conclude, we have found that the DM interpretation of the PAMELA positron excess can be tested by the FERMI 
LAT telescope in the diffuse emission at mid-latitudes and high energy. In the case of a DM candidate injecting a 
very hard e+/e~ yield, the signal is given by the IC -I- FSR emission. The case of a DM candidate with a sizable and 
hard 7-ray yield from 7r°-decay (like the case of r^/r" as final state of annihilation/decay), can be even more easily 
tested. A crucial ingredient for the discussion is the different spatial profile for the CR primary sources (confined 
to the disc) and of dark matter induced components (spherical distribution). The two terms can be disentangled 
looking at the angular profile of the diffuse emission, with the optimal region to single out the DM component being 
at intermediate latitudes. Such cross-correlation test of the PAMELA excess would then be performed by focusing 
on a nearby portion of the Galaxy, where the extrapolation on the DM density profile as well as on the propagation 
model parameters (from the locally-measured CR spectra) can be regarded as rather robust. Indeed, would FERMI 
discover an extra 7-ray term, with the spectral and angular features matching the features we discussed for the DM 
source, this would be a solid step towards the identification of the DM component of the Universe. 

In turn, would FERMI find that the 7-ray diffusion emission is instead in agreement with the prediction from 
standard CR components only, tight constraints on the DM interpretation of the PAMELA positron excess would 
follow. On the other hand, such picture would not be in contradiction with other scenarios addressing the PAMELA 
excess. If the additional e^/e~ sources have a spatial distribution analogous to the standard CR source distribution, 
such as, e.g., in the case of pulsars or secondary production inside CR sources, it will be very hard to single out the 
associated Galactic diffuse emission. Indeed, it would be, typically, subdominant and with a signal to background 
ratio at the same level in any portion of the sky; hence, a different strategy would be needed to cross-check these 
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interpretations of the PAMELA positron excess. 
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